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Abstract 


The objective of the present paper is to investigate the effect of clamping force on the performance of proton exchange membrane fuel cell 
(PEMFC) with interdigitated gas distributors considering the interfacial contact resistance, the non-uniform porosity distribution of the gas diffusion 
layer (GDL) and the GDL deformation. The clamping force between the GDL and the bipolar plate is one of the important factors to affect the 
performance and efficiency of the fuel cell system. It directly affects the structure deformation of the GDL and the interfacial contact electrical 
resistance, and in turn influences the reactant transport and Ohmic overpotential in the GDL. Finite element method and finite volume method are 
used, respectively, to study the elastic deformation of the GDL and the mass transport of the reactants and products. It is found that there exists an 
optimal clamping force to obtain the highest power density for the PEMFC with the interdigitated gas distributors. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


The proton exchange membrane fuel cell (PEMFC) system 
is considered as an alternative power source for electric vehicles 
by virtual of its high-energy efficiency, pollution-free charac- 
teristic, and simplicity in design and operation. It has received 
more and more attention in the world during the last decade. 
Numerical and experimental studies have been rapidly increased 
in order to improve the performance and cost-effectiveness. 
Experimental results have shown that a suitable design of the 
interdigitated flow field gives rise to a great improvement of the 
system performance [1]. The improvement is attributed to the 
fact that it converts the transport of the reactants and products 
in the porous electrode from a diffusion mechanism to a forced 
convection mechanism. Furthermore, the shear force exerted by 
the gas flow helps to remove most of the liquid water from the 
porous electrode, thereby minimizes the possibility of the elec- 
trode flooding. 

Several attempts of mathematical modeling have been made 
to understand the advantages of the interdigitated flow field 
design. Yi and Nguyen [2] presented a multi-component model 
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of the hydrodynamics of the reactant gas in the porous elec- 
trode of the PEMFC with the new design flow field. This model 
was then extended to the two dimensional and two-phase multi- 
component transport model by He et al. [3] to account for the 
effects of liquid water in the porous electrode on the PEMFC 
performance. They investigated the effects of the operating con- 
ditions and the design parameters of the air cathode on the 
performance of the fuel cell. The study indicates that the sat- 
uration of liquid water at the cathode greatly affects the system 
performance and efficiency. Using the similar model, Wang et al. 
[4,5] investigated effects of the magnetic particles on the perfor- 
mance of the PEMFC with the interdigitated gas distributors. All 
the models mentioned above have assumed for simplicity that 
the porosity of the GDL is constant, i.e., the GDL is assumed to 
be arigid material whose porosity is independent of the camping 
force. 

In order to seal the fuel cell against any gas leakages and 
reduce the interfacial contact resistance, the gas diffusion layers 
and the bipolar plates are usually bolted together under a suitable 
pressure, as shown in Fig. 1. Obviously, the pressure will reduce 
the porosity below the shoulder area, which in turn increases the 
convection resistance of the reactants and products. The exper- 
imental observation [7] has shown that the clamping force will 
strongly affect the power density of PEMFCs. Roshandel et al. 
[8] investigated the fuel cell performance considering the effects 
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Fig. 1. Schematic of a PEM fuel cell stack. 


of porosity variation in the GDL. Their numerical analysis show 
that a decrease in the porosity reduces the reactant transport in 
the GDL. Chu et al. [9] investigated effects of the non-uniform 
porosity caused by liquid water on the fuel cell performance in 
terms of some representative parameters, such as oxygen con- 
sumption, power density, etc. The non-uniformity of the porosity 
in the GDL is supposed to be described by four different con- 
tinuous functions of the coordinate position. However, all of 
those models considering the non-uniformity of the porosity 
have assumed that the porosity varies as a simple function given 
previously. Zhou et al. [6] recently studied the interfacial con- 
tact resistance and elastic deformation of the GDL using finite 
element method, and found that the porosity distribution of the 
GDL strongly depends on the clamping force and the structure 
of the bipolar plates. It is found that the porosity distribution 
exists in a complex manner, but cannot be described by a simple 
function given previously. 

The objective of the present work is to study the effects of 
the compression pressure due to the clamping force on the per- 
formance of PEMFC with the interdigitated gas distributors. 
In this work, we first use a finite element method (FEM) to 
analyze the contact pressure, the deformation of GDL, and the 
porosity distribution, then we analyze the mass transport using 
a two-phase flow model. The numerical results show that a high 
clamping force results in a lower porosity in the GDL, which will 
increase the transport resistances of both the gas phase and liquid 
phase. On the other hand, a high clamping force may decrease 
the interfacial contact resistance, and hence minimize the elec- 
trical resistance loss inside the fuel cell. These two effects may 
work together to yield an optimum clamping force for the fuel 
cell. 


2. Model description 


The computational domain, i.e., the region surrounded by the 
dashed line as shown in Fig. 2a, consists of the gas diffusion layer 
and a period of the interdigitated gas distributor. Although the 
studied domain of the FEM model to analyze the GDL deforma- 
tion can be half of what is adopted in the two-phase flow model 
considering the symmetrical characteristic of the GDL deforma- 
tion, we still use the same domain as the two-phase flow model 
in order to easily transfer the computational data between the 
two models. 


2.1. FEM model 


The FEM analysis model, as shown in Fig. 2b, is used to 
analyze the elastic deformation of the GDL and the porisity 
destribution. It is assumed that the cathod and the anode of 
the bipolar plates have the same structure, and the compressive 
modulus of the catalyst and membrane layer (about 100 MPa, 
T65-RH30) [10] are much higher than that of the gas diffu- 
sion layer (~10 MPa) [11]. Therefore, the deformation of the 
catalyst and membrane layer can be neglected, and will not be 
considered in the FEM analysis of the elastic deformation of 
the cell structure. Normal displacements at the left, right and 
bottom boundaries of the gas diffusion layer were restricted 
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Fig. 2. A schematic of (a) electrode of polymer electrolyte membrane (PEM) 
fuel cell and (b) FEM analysis model. 
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Table 1 
Parameters of the GDL used in the FEM analysis 


Porosity Ee (cross-plane) 
(MPa) 


E; (in-plane) 
(MPa) 


Interfacial resistivity, o 
(mQ cm?) 


0.7 10 100 p=2.0 (5.0p/E.)~°” for case A 


p=20 (2.0p/E-)~! for case B 


The interfacial resistivities for cases A and B are based on the experimental 
results using the graphite shoulder and the stainless steel shoulder, respectively 
[11,13,14]. 


due to the symmetry of the structure. Bipolar plate is treated 
as rigid structure. The vertical displacement of the shoulder is 
given downwards, simulating the stacking process. The problem 
mentioned above can be described as a large deformation contact 
problem with plane-strain boundary [12]. Plane-strain elements 
were used and contact elements were built at the possible contact 
region, as shown in Fig. 2. 

We use carbon cloth or carbon paper as the GDL material. 
The carbon paper is a kind of non-woven paper of carbon fibres, 
which is an anisotropic porous material. The mechanical prop- 
erty of the GDL and the relationship between the contact pres- 
sure and the interficial electric resistivity should be considered 
to obtain the porosity distribution and total interficial resistance. 
However, only a few open literatures have been reported on the 
modeling of the mechanical performance of the GDLs and the 
total contact resistance between the GDL and the shoulder. The 
parameters employed in the present paper are given in Table 1. 
Most of the data are based on the experimental results and simu- 
lating models of the similar material [11,13—15]. Cases A and B 
are based on the experimental results using the graphite shoulder 
and the stainless steel shoulder, respectively. The Ohmic overpo- 
tential depends on the parameters used in the contact resistance 
model (see Table 1). However, we found that they will not influ- 
ence the general conclusions if the model of contact resistivity 
does not change but the contact resistivity varies in the same 
order. 


2.2. Porosity distribution and interficial resistance 
calculation 


The GDL deformation and the interfacial contact pressure 
are first obtained using the FEM analysis. Then, the porosity 
distribution and interficial resistance are calculated. Based on the 
large deformation theory of solid mechanics [12], the porosity 
after stacking is: 


o= Vp — Voo E Vod =e") E 0 -1 +e®V 0) 


V Voetv eéV 


where Vo is the initial total volume of the porous media, ¢o the 
initial porosity, and £y the bulk strain, Vp is the pore volume. 
The top-bar indicates the parameters after stacking. The volume 
of the porous media is calculated as: 


P = Voe = Voe tte: D 


The total interfacial contact resistance in a period of the fuel 
cell is [6] 


-1 


1 n 
= 2 X Ai/ pi (3) 
i=l 


R=2 Sor 
i=1 


where R;, Aj, and p; are the resistance, contact area, resistivity of 
the ith contact element, respectively. This equation is based on 
the assumption that the contact surfaces are equipotential. There- 
fore, the interfacial contact resistance between the two surfaces 
is considered as an electric circuit in parallel connections. The 
factor 2 in Eq. (3) is due to the same contact resistances in the 
anode and cathode. 

The total compression force in the unit width (z direction), 
F, has the following relationship with the contact stress compo- 
nents: 


n 
F=) pint pli (4) 
i=] 


where pi, is the normal contact stress component (contact pres- 
sure), A the tangential contact stress component (friction), 7 
and T are the unit vectors in normal and tangential directions of 
contact elements, respectively. The nominal clamping pressure, 
Dc, is defined as: 


F 
=> (5) 


À 
where A is the length of the period of the cross-section of the 
bipolar plate as shown in Fig. 2a. 


2.3. Two-phase flow model of the cathode 


The following assumptions are used in the two-phase flow 
model of the cathode for studying the effect of the clamping 
force on the performance of PEMFCs: (1) steady flow state; (2) 
isothermal flow; (3) laminar flow in the GDL; (4) the gas and 
liquid in the electrode flow as continuous phases described by 
Darcy’s law; (5) the water is generated as liquid at the membrane 
and electrode interface by the oxygen reduction reaction and 
the net water migration from the anode. Based on the above 
assumptions, the continuity equations are: 


V -(C8v®) + ry =0 (6) 


for the gas phase, and 
—ry =0 (7) 


for the liquid phase, where C® is the gas molar concentration, vē 
the gas superficial velocity, v! the liquid superficial velocity, rw 
the interfacial mass-transfer rate of the water between the gas 
and liquid phases, py the water density, My the water molecular 
weight. 
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The momentum conservation are given in the following 
according to Darcy’s equation: 


KE 
v=- (5) VPE (8) 
ps 


for the gas phase, and 


K! 
v=- a VP (9) 
H 


for the liquid water, where P®, uë and K F are the pressure, 
viscosity, effective permeability of the gas, respectively. P!, u! 
and K a are the pressure, viscosity, efficient permeability of the 
water, respectively. 

The species conservation equations are: 


V: (Cov! - D&-VCo) =0 (10) 
V = (Cwo® — DV Cw) + rw = 0 (11) 


for O2 and H20 vapour, respectively. Where DÌ; is the effec- 
tive diffusivity of oxygen, DY, the effective diffusivity of water 
vapour, Co the molar concentration of oxygen, and Cw the molar 
concentration of water vapour. 

The interfacial mass-transfer rate of the water between the 
gas and liquid phases, rw, can be expressed as: 


rw = ke@(1 — s)yw( RT)! (yw P® — PS#) q 
+ kypspwM,,' (yw P® — PS") (1 — q) (12) 
with a switch function 


gx Lt iw PE = PHD PE = PV) 
2 


(13) 


where yw = Cw/C8, s is the saturation of the liquid water in the 
pore, T the temperature (°C), R the gas constant, Mw the molec- 
ular weight of water. When the water partial pressure (ywP®) 
is greater than the water saturation pressure (PS), the switch 
function q = 1. Otherwise, q =0. 

The boundary conditions have been described in the litera- 
tures [3,16], which simulate the PEMFC with the interdigitated 
flow field. At the boundary where x=0 or x=A (see Fig. 2b), 
the symmetry conditions are used for all the parameters, and the 
normal gradients of the four independent variables, P8, Co, Cw, 
and s, in governing Eqs. (6)—(13) vanish. The same boundary 
condition is also employed at the interface of the GDL and the 
shoulder of the gas distributor. At the GDL/catalyst-layer inter- 
face, itis assumed that there is no normal superficial gas velocity. 
The mass fluxes of oxygen and liquid water can be defined as 
the functions of the local current density Z as follows: 


aC I 
-Dp => 14) 
oy y=0 4F 
Ply 1 I 
E, SEON ee Z 15 
My” y=0 +a) F i ) 


The liquid water flux at the interface arises from two sources: 
the oxygen reduction and electro-osmosis. The parameter a, 
which is assumed to be constant, represents the net water trans- 
port coefficient in the membrane. The mass flux of the water 
vapour at the interface equals zero, i.e., 


pw aCw 


y=0 


The simplified Butler—Vollmer equation is used to describe 
the local current density along the catalyst-layer as a function 
of the oxygen concentration at the GDL/catalyst-layer interface 
and the activation overpotential 7 


@ Co (AF) 
I= l exp 
0 Co, ref RT 


where Jp is the exchange current density, Co,ref the oxygen 
reference concentration, k the transfer coefficient of oxygen 
reduction reaction. The two-phase flow model described above 
is similar to that presented by He et al. [3]. However, at the 
present work, some parameters and equations are described more 
accurately. 

In the governing equation, the effective permeability and dif- 
fusivity are both correlated with the porosity and the saturation 
level of the liquid water. In addition, the diffusion coefficient is 
a function of temperature and pressure. Tomadakis and Sotir- 
chos [17], and Nam and Kaviany [18] numerically simulated 
the effective diffusivity for a random fibrous porous media, and 
suggested that: 


Dest = Dh(T)f ($) g(s) 


_ T\". $ — bp 2 b3 
= (7) 6(f=e) a-s) a8) 


where ġp is a percolation threshold, D the diffusivity at T, when 
p= 1 and s=0, and by, b2 and b3 are empirical constants. bı is 
2.334 and 1.832 for vapour water and oxygen, respectively. For 
porous media composed of two-dimensional, overlapping, and 
random fibre layers: p =0.11, b2 is 0.521 and 0.785, b3 is 2 
and 4, for in-plane and cross-plane diffusion, respectively. The 
values of those parameters are based on the numerical analysis 
results considering the randomly prescribed pore-filling water 
distribution, and can be used in the analysis of fuel cell system 
with a carbon fibre diffusion layer [16,17]. 

In the present work, the well-known Carman—Kozeny equa- 
tion supported by experiments [19,20] is applied to describe the 
relationship of the porosity and the permeability for an individ- 
ual phase flow. We assume that the permeability of the gas phase 
gives the following relation with the phase saturation: 


(17) 


— (1 - 8 (19) 
1-¢ 


(1 y ~~ o 
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for cross-plane permeability. The assumption of the relationship 
of the efficient gas permeability and liquid water saturation is 
based on the research work of Chen et al. [21] and Schulz et al. 
[22], and the effect of liquid water saturation on the efficient gas 
diffusivity [16]. 

The permeability of the liquid water phase is usually 
described by a power function of its saturation, with an exponent 
ranging from 1 to 3 [23]. Most of the literatures use 1.5 and 3 in 
the simulation of PEMFC. The present paper takes it as 1.5, and 
gives the following expression for the permeability of the liquid 
water phase: 


(21) 


Although the Carman—Kozeny equation has some limitations 
to predict the permeability of carbon papers, many researchers 
show that Carman—Kozeny equation is still applicable if the 
Carman—Kozeny constant, K}, is redefined or treated as an 
empirical constant based on the experimental results [19,24,25]. 
Rodriguez’s experimental observations [19] showed that the 
Carman—Kozeny equation can well describe the permeability 
of the fibre mats (similar to a carbon paper) when the poros- 
ity ranges from 0.5 to 0.8. At the present work, the porosity 
ranges from 0.54 to 0.7, giving a reasonable range to use the 
Carman—Kozeny equation. 


2.4. Unsaturated flow theory (UFT) 


Unsaturated flow theory has been widely used in the modeling 
of fuel cells [3,26]. The pressure in the liquid phase is given by: 
P! = PE — pe (22) 


where P* is the capillary pressure. Eqs. (9) and (22) are combined 
to yield: 


1 l 
v= -Eor VPE + Ketty pe 
m H 
Kle ë KS K} OPS 
e a ype sot we (23) 
ul Ke. us ul as 


The capillary pressure between the gas and liquid phases is 
related to the saturation of liquid phase [27]: 


1/2 
a ) Hs) (24) 


PE =o cos(@&) ($ 


where J(s) is the Leverett function, ø the surface tension for 
liquid water and air system, 6, the contact angle of the GDL 
surface. For a hydrophobic GDL [28], 


J(s) = 1.4175 — 2.120s* + 1.2635? (25) 
Making use of Eqs. (19)-(21), and (23)-(25), one obtains: 


1/2 
1_ Khe sto cos(@)(Kh) $ as 
yl KE (=s) do 1-9 
x (1.417 — 4.24s + 3.7882)Vs (26) 


where n is 2 and 4 for in-plane and cross-plane permeability, 
respectively, which has been introduced in Eqs. (19) and (20). 
In this study, the gas and liquid phases within the GDL exist as 
continuous phases, as a result that the Darcy’ law is used to solve 
the transport phenomena of gas and liquid water, respectively. 
Eqs. (21) and (23)-(25) come from the literatures of Wang and 
co-workers [26—28] and Nguyen and co-workers [3]. 


2.5. Calculation of fuel cell voltage 


The polarization curve of a PEM fuel cell plays an important 
role in determining the performance of the fuel cell system. Due 
to the irreversible energy loss in a practical fuel cell, the actual 
potential (V) is decreased from its open circuit voltage (E). These 
losses include activation overpotential (AVact), Ohmic over- 
potential (AVohm), and concentration overpotential (A Vconc). 
Therefore, the operational voltage is [29]: 


V=E-— AVohm E A Vact = AVeone (27) 


The open circuit voltage, E, is decreased from its thermody- 
namic potential due to gas crossover. The Ohmic overpotential 
due to the resistance of the membrane and ions in the cell com- 
ponents (only the membrane resistance and contact resistance 
are considered here) is given by [29]: 


AVohm = Z (Rcr T Rmemb) (28) 


where Rep is the contact resistance between the GDL and the 
bipolar plate, and Rmemb is the flow resistance of H* ions in 
the electrolyte membrane. The activation overpotential is given 
by the Tafel equation, which is a transformation of Eq. (17). 
Concentration overpotential is given as a function of the limiting 
current in the following form [29]: 


AV, - lnį| 1 Í 29 
cone = 7 n( =) (29) 
where Jim is the limiting current density. In reality, however, 
the assumption of a completely mass-transfer limiting case is 
rarely valid because there is a concentration dependence in the 
activation kinetics of reaction that affects activation polarization 
as well [30]. Thus it is neglected at the present simulation. 


2.6. Solution strategy 


The governing equations of two-phase flow model described 
in Section 2.3 are discretized using a finite volume staggered 
scheme. The code is designed using Fortran Language. The well- 
known SIMPLE algorithm and upwind scheme are utilized in the 
numerical solution. The studied domain is divided into rectangle 
elements with four nodes. Finite element analysis for the elastic 
deformation of the GDL is first carried out to obtain the real 
shape of the GDL after deformation. The new grid is changed 
into the approximate orthogonal curvilinear grid. This is helpful 
to the discretization of the governing equations as described in 
Section 2.3. 

The effect of the number of elements on the computational 
results is investigated. Numerical tests show that 1600 (20 in 
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the height and 80 in the length) elements provide a sufficient 
solution accuracy, giving a good agreement with the numerical 
result of the rigid GDL [16]. Our numerical solution gives a 
very good numerical stability, but never meet a convergence 
problem. 


3. Numerical results and discussions 


The operating and design parameters used in the present 
simulation are listed in Table 2. The distributions of the oxygen 
mole concentration and the liquid water saturation with different 
clamping pressures are shown in Fig. 3. Effect of the clamping 
pressure on the performance of the fuel cell is considerable large 
based on our numerical analysis. To the best of our knowledge, 
however, no literatures are reported before to consider the 
deformation of gas diffusion layer caused by the clamping force. 
The porosity and thickness of the GDL are decreased due to the 
compression generated by the clamping force. The transport 
resistance increases with decreasing diffusion and permeability 
coefficients depending on the GDL porosity. Fig. 3a and Fig. 3b 
display the distributions of the oxygen mole fraction and the 
saturation level without considering the elastic deformation of 
the GDL, but Fig. 3c and Fig. 3d show the different distributions 
when the elastic deformation of GDL is considered. Obviously, 
there is a high saturation of liquid water near the interface of 
the catalyst-layer and the GDL below the inlet due to the GDL 
deformation (see Fig. 3d). This high saturation is caused by the 
low transport speed of the gas phase. The friction drag force on 
the interface of the gas and liquid is helpful to the removal of the 
liquid water. The saturation of the liquid water leads to a small 
effective reaction area. The local current density distributions 
corresponding to Fig. 3 are shown in Fig. 4 for an elastic GDL 
and a rigid GDL, respectively, when the nominal clamping 


6.036 


4 842 


879 
Table 2 
The basic parameters used in the numerical analysis 
Parameter Value 
Inlet channel width (half)* 0.50mm 
Shoulder width* 1.00 mm 
Round radius of shoulder angle? 0.05 mm 
Outlet channel width (half)* 0.50 mm 
Electrode height* 0.25 mm 
Membrane thickness (Nafion 112)° 0.051 mm 
Gas permeability of the electrode, KP 1.0 x 1071? m? 
Liquid-water permeability of the electrode, K G 1.0 x 107!? m? 
Dry porosity of the electrode, po? 0.7 
Inlet mol fraction of oxygen 0.21 
Inlet mol fraction of water 0.00 
Inlet mol fraction of nitrogen 0.79 
Inlet channel pressure? 1.007 atm 
Outlet channel pressure? 1.0 atm 
Temperature* 60°C 
Gas viscosity, 2° 2.059 x 1075 Pas (60°C) 
D? (at T,)* 1.775 x 1075 m? s7! (0°C) 
DY (at Ty)? 2.56 x 1075 m? s7! (34°C) 


Condensation rate constant, k,* 
Vapourization rate constant, ky* 

Net water transport coefficient, œ? 

o cos(8e)? 

Liquid-water viscosity, u! 

Proton conductivity of membrane, om° 


100s7! 


100 atm! s7! 
0.5 
—0.0001 Nm7! 


4.7 x 1074 Pas (60°C) 
1027! m7! 


a Reference [3]. 
b Assumed. 
© Reference [31]. 


pressure is 1.93 MPa (the maximum normal deformation of 
GDL is 0.0753 mm). The average current density considering 
the GDL elastic deformation is less than half of the current 
density without considering the GDL deformation (rigid DGL). 
This is the reason why we should pay more attention to the 
GDL deformation caused by the clamping force. In fact, the 


1.261 


MIN 


MIN LO.068 


4 792: 3.6013] 


(d) 


Fig. 3. Comparison of oxygen and water distribution of the fuel cell: (a) the oxygen mole concentration and (b) the saturation level of the liquid water without 
considering the GDL deformation; (c) the oxygen mole concentration and (d) the saturation level of the liquid water considering the GDL deformation. 
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Fig. 4. Current density distribution along the catalyst-layer/GDL interface in 
the cathode-length direction for the elastic gas diffusion layer (dashed line) and 
the rigid gas diffusion layer (solid line). The nominal clamping compression 
Pc = 1.93 MPa. 


GDL deformation leads to the change of GDL thickness, except 
for the change of porosity distribution. The performance of 
the fuel cell with the interdigitated flow field is sensitive to 
the GDL thickness [3]. Grujicic and Chittajallu [16] made an 
optimization for the PEMFC cathode geometrical parameters 
and concluded that there is an optimum thickness of GDL. 

The polarization curves considering different potential losses 
are plotted in Fig. 5. In the region of Ohmic polarization, the 
potentials due to the different contact resistances for the graphite 
and stainless steel shoulders are completely different. When the 
nominal clamping pressure is 1.93 MPa, the Ohmic resistance 
of the membrane, the contact resistances with the graphite and 
stainless steel shoulders are 2.54, 0.252 and 5.2 mQ m, respec- 
tively. The clamping force is an important parameter to control 
the Ohmic loss for the shoulder having the contact resistivity 
with the carbon paper as described in case B in Table 1, for 
example, the stainless steel bipolar plate. However, for the other 
types of shoulders having a small contact resistance with the 
carbon paper as described in case A in Table 1, for example the 
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Fig. 5. Polarization curves of the two-phase flow prediction with different shoul- 
der materials of the bipolar plates. 
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Fig. 6. Polarization curves for different nominal clamping pressures. 


graphite bipolar plate, too large a clamping force is not expected 
due to the higher transport resistance. 

Comparisons of the potential and power density of the fuel 
cell are given in Figs. 6 and 7 for the several clamping pressures 
as indicated. Effects of the nominal clamping pressure on the 
interfacial contact resistance and porosity are in two opposite 
ways: a high clamping force can give a low interfacial contact 
resistance as expected, but also gives rise to a low porosity of the 
GDL, which may decrease the system efficiency. For a high con- 
tact resistivity shoulder, the contact resistance plays a dominant 
role due to its large contact resistance at a low clamping force. 
But at a high clamping force, the contact resistance of this kind 
of shoulder and carbon GDL decreases slightly with the clamp- 
ing force. This time the GDL deformation plays a dominant role 
in determining the power density. However, when the shoulder 
with a low contact resistivity is used, the GDL deformation is 
an important factor due to the small contact resistance. Based 
on the analysis above, we can see that there must exit an opti- 
mum clamping force for a practical fuel cell system. The power 
density decreases with increasing the clamping force when the 
influence of Ohmic contact resistance is smaller than that of the 
GDL porosity. 
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Fig. 7. Power density vs. current density for different nominal clamping pres- 
sures. 
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At the small clamping pressures, what are shown in 
Figs. 6 and 7 are quantitatively in good agreement with the 
experimental measurements [32] for a micro-PEMFC with the 
interdigitated flow field. Unfortunately, at the high clamping 
pressures, we do not have any experimental data available to 
compare for the PEMFC with the interdigitated flow field. But 
good qualitative agreements occur at the polarization and power 
density curves predicted by the present paper for the interdig- 
itated flow field and by Chu et al. [9] for the one-dimensional 
model if the porosity change is considered, especially at the 
high current region. We can estimate that the clamping pres- 
sure affects the performance of the PEMFC in a similar manner 
for any types of flow fields. For the interdigitated flow field, the 
effect may be larger than the other types of flow fields based on its 
working mechanism. But this needs a further study. Lee et al. [7] 
measured the performances of PEMFCs with a serpentine flow 
field for three types of GDLs (ELAT®, TORAY™, a combina- 
tion of TORAY™ and CARBEL®) and found that the optimal 
clamping force is different for different GDLs. However, the 
changing fashions of the polarization and power density curves 
are in good agreement with what is predicted in the present 
paper. 

The current density corresponding the maximum power den- 
sity in the power density-current curve depends on the slope 
of the polarization curve in the Ohmic region, as shown in 
Figs. 6 and 7. The slope is controlled by the Ohmic loss and 
the mass-transport resistance which are affected, respectively, 
by the contact resistance and the increased transport resistance 
of reactant caused by the clamping force. 


4. Conclusions 


Finite element method (FEM) and two-phase flow model 
are used to investigate the effect of the GDL deformation 
caused by the clamping force on the performance of an air 
cathode of the PEM fuel cell with the interdigitated gas distrib- 
utors. From the numerical results, the following conclusions are 
reached: 


(1) The clamping force affects the permeability and diffusion 
of the reactant gas and the transport of the liquid water due 
to the GDL deformation and the porosity variation. A low 
clamping force increases the interfacial electrical resistance, 
but a high clamping force decreases the GDL porosity and 
thus increases the transport resistance of the gas and liquid 
water. 

(2) For a practical fuel cell system, there would be a maximum 
power density if an optimal clamping force is found. The 
optimal design parameters for areal fuel cell system depends 
on the mechanical and electrical properties of the bipolar 


plates and gas diffusion layer, and the interfacial electrical 
property and the across shape of the bipolar plates, etc. 
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